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ABSTRACT 


A systematic analytical approach to the determination of helicopter 
IFR precision approach requirements Is formulated. The approach Is based 
upon the hypothesis that pilot acceptance level or opinion rating of a 
given system Is Inversely related to the degree of pilot Involvement In 
the control task. A nonlinear simulation of the helicopter approach to 
landing task Incorporating appropriate models for the UH-1H aircraft, 
the environmental disturbances and the human pilot was developed as a 
tool for evaluating the pilot acceptance hypothesis. The simulated pilot 
model Is generic In nature and Includes analytical representation of the 
human Information acquisition, processing, and control strategies. Slmul 
tlon analyses In the flight director mode Indicate that the pilot model 
used Is reasonable. Results of the simulation are used to Identify candl 
date pilot workload metrics and to test the well known performance-work- 
load relationship. A pilot acceptance analytical methodology Is formula- 
ted as a basis for further Investigation, development and validation. 
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INTRODUCTION 


Background 

It Is safe to state that helicopter operations are expected to grow 
rapidly over the next two decades [1). This optimistic forecast for ex- 
panded potential markets Is predicated upon the assumption of an established 
capability for all-weather rotorcraft operations, and the ensuing advant- 
ages of reduced direct operating costs. However, at the present time, 
helicopters under Instrument Meterologlcal conditions (INC) are forced 
to operate as though they were a fixed-wing aircraft and must abide by the 
standard Instrument Flight Rules ( IFR) regulations and weather minima. 

This approach does not give any credit to the helicopters most unique and 
distinguishing characteristic - Its ability to make steep decelerating 
approaches to a hover and touchdown Into confined landing sites. The 
pilot workload demand for a decelerated approach to a hover Is extremely 
high or not acceptable because of the generally poor handling qualities of 
the basic helicopter In low speed flight. Consequently, some minimum 
level of control and display sophistication In the form of a stability and 
control augmentation system (SCAS) and/or flight director guidance is, essen- 
tial for acceptable IMC precision decelerating approaches. 

There are a number of factors which Impact upon the pilot acceptance 
of a helicopter decelerating steep approach and landing (DSAL) task under 
IMC. Clearly, the most Important, though obvious, factor is the severity 
of the weather conditions as prescribed by the slant visibility and cloud 
celling at the landing site. Typically, weather minimums are given In 
terms of the runway visual range (RVR) and decision height (DH) associated 
with the standard ILS approach trajectory. However, the helicopter is 
uniquely qualified to perform steep precision approaches with glide slope 
angles as high as 12 degrees. Hence, more realistic approach weather 
minimums, specifically tailored to helicopter IFR approaches to a hover 
over a prescribed landing area (as opposed to a runway), need to be defined. 
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Depending upon the flight path glideslope angle Y re f and weather minima, 
some degree of deceleration under IMC conditions may prove necessary. For 
example, helicopter stralght-ln approaches to a helipad under IMC to a 
DH of 50 ft would provide a horizontal distance to the landing site of. ap- 
proximately, 1000 ft for a 3* approach (DH/tan Y re f for a y approach 
and decision height DH) In which to decelerate from breakout (IFR/VFR 
transition) speed to hover. Results of a flight Investigation of visual 
approaches to a hover, conducted by NASA-LaRc [2], show that the pilot 
chooses to decelerate In two phases* 1) a slow deceleration, starting at 
a range of 3000 ft, from an Initial groundspeed of 50, 80 or 100 knots 
to, approximately, 45 knots at a range of 1000 ft, followed by 2) a more 
rapid deceleration to a hover on the prescribed landing pad. No such 
data exists for the helicopter IFR approach to a hover task. However, 
the VFR flight test data does provide some Indication of the close rela- 
tionship between aircraft groundspeed and dlstance-to-cjo from the landing 
pad. Thus, the need for deceleration under IMC conditions would depend 
upon the glideslope angle of the approach profile and the prescribed de- 
cision height for IFR/VFR transition. Pilot acceptance would be contin- 
gent upon his being able to perform the landing task with satisfactory 
performance and within reasonable workload demands. 

The objective of helicopter IFR operations Is to try and reduce the 
weather mlnlmums to as low a level as possible (zero vlsiblllty/zero 
celling). For given weather mlnlmums and aircrew manning level, the 
acceptability of a helicopter for IFR approach to a landing pad would 
depend upon the quality or sophistication of the guidance, navigation, 
control, and display components of the overall system. At the present 
time, the acceptance of a given helicopter with its inherent dynamics, 
flight control system and display avionics Is based upon subjective 
flight evaluation by the examining pilot on an Individual basis. A pilot/ 
aircraft system Is accepted for a given IFR approach scenario If It can 
meet prescribed performance standards while keeping crew pilot workload 
levels within reasonable bounds. Such an approach Is Intuitive In nature 
and provides little Insight Into the rationale for the examining pilots 
acceptance criteria. Therefore, there Is a clear need to develop a 
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rigorous structured approach towards the evaluation and, If possible, 
specification of aircraft system requirements (l.e., crew, displays, 
flight control, navigation and guidance) for helicopter IFR precision 
approach and landing operations. 


Objectives 

The NASA Ames Research Center has Initiated a long-term program 
aimed towards the development of a technology base tor Improved heli- 
copter all-weather capability. A step In pursuit of this goal was to 
conduct systematic helicopter piloted simulation experiments for a pre- 
cision decelerating approach to a hover for various levels of flight- 
control and display sophistication. The goal of these experiments 
(referred to as the single-pilot IFR or SPIFR simulation) was to provide 
an adequate data base for: 


1. the experimental evaluation of minimum guidance, control, and 
display requirements for various weather minima; and 

2. the formulation of a systematic analytical methodology directed 
towards the establishment of minimum guidance, control and 
display requirements for acceptable, single-pilot, helicopter 
precision landing approaches, under IMC. 

A description of the SPIFR simulation Including the analysis of the 
experimental data Is given In reference 3. The objective of this research 
effort Is to use data acquired In the previous fixed base simulation In- 
vestigation of helicopter SPIFR requirements towards the formulation of 
an analytical methodology for determining pilot acceptance. The overall 
approach Is discussed In the following paragraphs. 


Technical Approach 

The current state-of-the-art for assessing competing aircraft system 
configurations Is fundamentally empirical. The level of acceptance or 
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satisfaction Is usually expressed In terms of a pl.ot opinion rating (POR) 
scale [4,5] such as the Cooper-Harper handling qualities scale. For lack 
of a better alternative the POR approach Is the most frequently used method 
for evaluating the relative merits of candidate aircraft system configura- 
tions. 

The degree of pilot acceptance expressed as a pilot opinion rating 
depends upon a number of factors which can broadly be put Into two cate- 
gories • system performance and pilot workload. In axiomatic terms, a 
system A would be rated better than system B, If for a given acceptable 
level of system performance, system A requires lesr. pilot workload than 
system B. Clearly, a tradeoff between the two attributes of system per- 
formance and pilot workload is Involved In arriving at an aggregate 
measure such as an opinion mating of pilot acceptability. 

Further quantification requires the development of analytical measures 
for system performance and pilot workload. This constitutes the funda- 
mental step and stumbling block In the formulation of an analytical methodo- 
logy for determining the level of pilot acceptance or opinion rating. 

The problem Is complicated further by virtue of the fact that In most real- 
istic conditions, system performance Is a vector quantity and that pilot 
workload may also depend upon many Individual factors. However, the system 
performance vector 

s - (Sj. s 2 , . . , s n ) T (1) 

Is relatively easier to quantify In terms of objective or physical attri- 
butes, given adequate knowledge and background of the particular phase 
of the piloting task. Thus, for a helicopter precision hover task, the 
system performance vector may be selected to be 

s - (o x , o q , o y , o p , o z , o*) (2) 

where o^ represents the root mean squared (rms) deviation of variables 
(x*q,y,p,z,z) from nominal, namely (x o ,0,y o ,0,z o ,0). This representa- 
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tlon for the system performance vector ts not u.ilque, and could Include 
standard deviation of additional variables such as sideslip 6 and yaw 
rate r. Depending upon the task, the pilot must somehow combine the In- 
dividual components of the system performance vector to formulate an 
aggregate scalar metric, J, for total system performance; thus 

0 ■ 9(s) 

The functional g Is Internal to the human pilot and must be Inferred using 
systematic methods: for example - pilot questionnaire, multiattribute 
utility analysis, stepwise regression analysis and others. Typically, 
the functional g Is chosen In an ad-hoc manner based upon subjective 
judgement and pilot comments, and has the additive structure: 

J ■ OjSj ♦ OJ S 2 ♦ . . ♦ a„S n 

with the {<Xj} specified by the analyst. Nevertheless, Inspite of the 
arbitrariness of this procedure. It Is conceivable In terms of physical 
attributes such as standard deviations of key variables and hence amen- 
able to systematic analysis. 

On the other hand, finding the elements of the vector workload 
metric and the manner In which they contribute to an aggregate measure 
or Index for pilot workload Is even less obvious. The term workload has 
various connotations In the human factors literature. Typically workload 
has been partitioned into three categories: (1) scanning or information 

acquisition workload, (2) central Information processing or decision work- 
load, and (3) motor or control workload. Three approaches have been pre- 
sented in the literature for quantifying pilot-workload. They are (1) 
subjective opinion ratings, (2) reserve attention capacity as measured by 
secondary task performance, for fixed main task performance, and (3) 
physiological indicators such as galvanic skin resistance, pupil dilation, 
heart-rate, electroencephalogram (eeg) features and others. However, 
these techniques are not useful as analytical measures of workload. As 
a ■ onsequuence, techniques based on an analytical representation or model 
of the human pilot's Information acquisition, processing and control 
strategy have been sought over the past three decades. 


Model -based measures by definition assume that an input-output mathema- 
tical model for the pilots control strategy exists and can be determined. 
Given that a pilot model can be found, the following axioms are used Im- 
plicitly In the formulation of workload metrics: 

1) For a given task, It Is hypothesized that there exists a nominal 
pilot strategy (for information acquisition, processing, and 
utilization) characterized by a nominal set of parameters which 
corresponds to minimum workload for a given acceptable level of 
performance, 

2) Workload Increases for deviation of the pi lot's strategy from 
nominal or, In other words, with the deviation of the model para- 
meter values from nominal, and 

3) Workload is related to the sensitivity of the system performance 
to perturbations of the model parameters from nominal values. 

Thus, analogous to the earlier discussion on system performance, pilot 
workload can be represented as a vector quantity 

0 * (0j. e 2 . . . , e m ) T (5) 

where the components { 0 ^ } represent the parameters corresponding to the par- 
ticular pilot strategy model parameterization. The aggregate scalar work- 
load metric W is assumed to be some function of 0; thus, 

W * h(e) (6) 

where h(e) could be some linear or nonlinear function on 6. 

Assuming that suitable scalar metrics for system performance J and 
pilot workload W have been developed, the pilot-opinion rating, R, can be 
formulated as 

R - f(J.W) (7) 
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As discussed earlier, some sort of tradeoff between performance J and work- 
load U Is essential In evaluating the relative merits of competing system 
configurations and arriving at a preference In terms of an opinion rating R. 
A plausible relationship between J, W and R Is Illustrated In Figure 1. 
where Iso-rating curves are drawn on a J versus W plot. The figure shows 
three constant rating curves Rj^, and R 3 corresponding to low (e.g., 3), 
medium (e.g., 5) and high (e.g., 7) pilot opinion ratings, respectively. 
Pilot workload Is bounded by W fl|ax on the high side and "min representing 
minimum value or open loop control. Similarly the system performance has 
a minimum value J m1n corresponding to maximum pilot workload levels and 
a maximum value that culminates with Instability of the pilot aircraft 
closed-loop system and subsequent loss of control. 



Figure 1. Performance - Workload Iso-Rating Contours - 
Illustrative Example 
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Thus 4 regions of pilot workload from low to extreme are obtained. 

For a given level of acceptable performance J > J Q , the horizontal line 
OABC defines different levels of pilot acceptance or option rating. Thus, 
If the three rating curves Rj.Rg, and R^ correspond to three different 
aircraft configurations under evaluation then system 1 would be preferred 
over system 2 which In turn Is chosen over system 3. Points A, B and C 
define the rating or preference level In terms of the workload levels 
and Wq required to maintain the prescribed level of performance 
J Q . Figure 2 shows the plot of ratings versus pilot workload metric 
W for various constant levels of performance J. 



Figure 2. Pilot Rating versus Pilot Workload 
Iso-Performance Curves 
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For a given acceptable level of system performance J, the pilot opinion 
ratings may thus be obtained from Figure 2 as a function of the workload. 


A generic description of a systematic methodology for determining 
pilot acceptance level or rating for competing system configurations has 
been presented. Formulation of the above pilot acceptance analytical 
methodology Is based upon one fundamental assumption - the existence of 
a suitable mathematical model describing the human pilot's Information 
processing and control strategy and an acceptable parameterization of 
the system performance and pilot workload metrics. 


Task Outline 


The purpose of this research effort Is to formulate an appropriate 
dynamic model for the pilot In the helicopter precision approach to hover 
task and evaluate the suitability of candidate model-based system perfor- 
mance and workload metrics using available SPIFR data. The ultimate 
goal is to use the analytical measures for system performance and pilot 
workload in formulating an assessment methodology for determining the 
level of pilot acceptance for any given system configuration. 

The helicopter IFR approach to landing task is a complex piloting 
task requiring nonlinear and time varying control laws on behalf of the 
human pilot. Thus, the use of two existing UH-1H approach and landing 
simulation programs was stipulated at the beginning of this study. The 
two programs are: (1) a UH-1H nonlinear simulation of the approach and 
landing task, and (2) a generic linear time varying covariance propaga- 
tion program useful for the propagation of helicopter approach path dis- 
persions during the precision landing task. 

The overall approach consists of four tasks, as follows: 
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Task 1. Automatic Mode Validation: 

The objective of this task Is to validate the approach performance 
predictions of the existing nonlinear simulation and linear time varying 
covariance propagation programs by comparison with ensemble statistics 
data for the automatic mode acquired from the SPIFR simulation. 

Task 2. Manual Mode Validation: 

The purpose of this task Is to validate the two existing programs 
using SPIFR data for the manual flight director mode of operation. Re- 
presentative pilot model structures are to be used in the flight director 
mode validation. Configurations involving raw data are outside the scope 
of this preliminary effort. 

Task 3. Evaluate Candidate Pilot Model Based Workload Metrics: 

The purpose of this task is to identify via simulation the existence 
of any plausible metrics for pilot workload and their relationship to 
system performance and pilot opinion rating. 

Task 4. Formulate Pilot Acceptance Analytical Methodology: 

Here the experienced gained from Tasks 1-3 above is to be assimilated 
to propose a systematic analytical approach towards determining pilot 
acceptance levels or opinion ratings. 

The above task delineation served the general purpose of providing 
a framework for this study effort. The approach followed this outline 
with some minor modifications dictated by the problems encountered during 
the course of this investigation. 
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PILOTED SIMULATION EXPERIMENTS 


I 

I 

I Facility Description 


The SPIFR experiments were conducted on the NASA-Ames V/STOLAND fixed 
base UH-1H simulator, shown In Figure 3, and Involved the participation 
of three pilots (2 NASA and 1 Am\y). Details on the V/STOLAND digital 
avionics system may be obtained from reference 6. A six degree of freedom 
mathematical model of a UH-1H helicopter [7] at a gross weight of 6158 lbs 
was used to drive the real time simulation. 


Figure 4 shows a schematic of the UH-1H simulator cockpit display 
panel. The Instruments are arranged in a standard "T" format except for 
the CRT multifunction display which was not used in the simulation. In 
addition, engine power Indicator Instruments (torque-meter and dual tacho- 
meter) and the inclinometer (turn and slip Indicator) were not operational 
during the experiments. 


*" ■ ■ ■ > 
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1. Attitude Director 
Indicator (ADI) 

2. Horizontal Situation 
Indicator (HSI) 

3. Barometric Altimeter 

4. Radio Altimeter 

5. Mode Select Panel 

6. Multifunction 
Display (MFD) 


7. True Airspeed 

8. Barometric Altitude 

9. Vertical Speed 
Indicator (IVSI) 

10. Torquemeter 

11. Dual Tachometer 

12. Clock 


FLR: Flare Annunicator - Green 

DH : Decision Height Annunicator - Amber 


Figure 4. V/STOLAND Cockpit Displays 


11 






ii? 



U 





The helicopter controls consisted of cyclic stick, collective stick 
and directional pedals. The controls travels are limited to + 6.3 In for 
the longitudinal cyclic, + 6.3 In for the lateral cyclic, + 3.45 In for 
the pedals and 11 In for the collective lever. Adjustable friction on the 
collective lever and artificial force -feel characteristics for the cyclic 
stick and directional pedals are also provided. 


Piloting Task 

The pilots were Instructed to follow the commanded flight path shown 
In Figure 5. The piloting task consisted of the following six phases: 

1) localizer capture, 2) tracking at cruise altitude h Q and constant speed 
V Q , 3) constant speed glide slope capture, 4) constant speed glide slope 
tracking until Initiation of deceleration or flare at X FLR* 5) IFR decelera- 
tion to velocity Vp H> at decision height, DH, and 6) VFR deceleration 
following breakout at DH to a hover 10 ft above the prescribed 125 ft dia- 
meter landing pad. 



Figure 5. Flight Path Geometry 


Constant attitude deceleration guidance (3,8] was used to compute the 
commanded velocity profile during the flare maneuver. Decision heights 
(OH ) of 100 ft, 50 ft and 0 ft were flown. The aircraft velocity at deci- 
sion height, Vq H was adjusted downward with decreasing decision height to 
45 knots (76 ft/s), 31 knots (52 ft/s), and 0 knots corresponding to the 
three decision heights of 100 ft, 50 ft and 0 ft respectively. This down- 
ward adjustment of with reduced DH was necessary to allow for the 
shorter distance available for an acceptable visual deceleration profile. 

A simplified Dryden wind model in each of the three aircraft body axes 
was used. A stationary model was selected through presimulation pilot 
comments since It is more amenable to analytical investigation as required 
In this study, and because of the pilot conments Indicating no noticable 
differences between a non-stationary wind model and the chosen stationary 
representation. 


Control and Display Configurations 

Four control augmentation levels were simulated as follows: 

1. Basic UH-1H with stabilizer bar (WSB) 

2. UH-1H (WSB) with control stick steering (CSS) 

3. UH-1H (WSB) with CSS and automatic collective axis 

4. UH-1H (WSB) with CSS and fully automatic control system. 

The Control Stick Steering (CSS) mode incorporates a first order 
attitude SCAS in pitch and roll axis, first order vertical velocity command 
control system (i.e., first order 6/ec, 4>/4>c and fi/fic systems), and a fully 
automatic directional (or pedals) control system. The CSS mode with auto- 
matic collective leaves only the pitch and roll axes under the pilot's 
control. The 4 axes, fully automatic system commands the aircraft to fly 
down the prescribed trajectory of Figure 5 to a hover following a constant 
attitude deceleration profile (CADP). 
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Two display configurations were simulated: (1) a raw data (RD) mode 
where the pilot had to use the standard cockpit Instruments In Figure 4 
and (2) a 3-cue flight director (FO) mode for the pitch, roll and vertical 
axes. The pitch and roll flight director signals drive the standard hori- 
zontal and vertical bars on the attitude director Indicator (ADI) shown 
In Figure 6, and are configured to provide a rate convnand/attltude hold 
system (i.e., k/s) between the pilots stick output and corresponding fllqht 
director display. The vertical flight director signal Is displayed as 
a "fly-to" symbol (except for pilot number 2 who preferred an error command) 
on the left side of the ADI. Two discrete signals on the ADI (see Fig. 6) 
are provided to aid the pilot - a green light to Indicate the Initiation 
of deceleration or flare and an amber signal to display the occurence of 
IFR/VFR transition at decision height. In the raw data display mode the 
pilot must use the standard glide slope and localizer Information for 
flight path tracking; IFR deceleration command Is missing and the pilot 
must use the raw data on altitude (radio altimeter) and airspeed on cue 
from the green light to arrive at the desired velocity at decision height. 
Note that the absence of a ground speed or range rate indicator on the 
cockpit panel forces the pilot to use airspeed and subjective estimates 
of the wind turbulence for Internally generating the deceleration velocity 
command. 


Experimental Test Matrix 

A total of 6 configurations corresponding to three levels of control 
augmentation and two types of displays were simulated. A fully automated 
system was also simulated to serve as a benchmark for evaluating competing 
piloted configurations and for validating the off-line UH-1H simulation 
programs 

Details of the specific system configurations tested are given in 
Table 1. Approximately 12 runs per configuration were conducted to provide 
an adequate data base for statistical evaluation of the results and 
subsequent use in the formulation of a pilot acceptance analytical metho- 
dology. 
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Data Summary 


Four types of data Mere obtained: (1) strip-chart recordings, (2) 
line-printer data, (3) digital tape recordings (run dum), and (4) pilot 
opinion ratings Including responses to a questionnaire. 40 variables 
were recorded on digital tape at a sampling period of 55 ms. Some of 
the key variables Include: aircraft attitude and angular velocities, 
airspeed, groundspeed, altitude, sink rate, cross-track error and error- 
rate, aircraft control positions, flight director commands, distance to 
helipad and longitudinal and normal accelerations. 

Cooper-harper ratings for the IFR portion alone, VFR portion alone 
and the entire IFR/VFR task were obtained. In addition, pilot ratings 
of the workload demands were also solicited. Details on the rating 
scales used may be found In reference 3. 

The raw data for the Individual runs were used to generate ensemble 
statistics (using time or distance to go as the Independent variable) for 
24 of the 40 recorded vralablts. These ensemble mean and standard devia- 
tion data provide the data base for validating the off-line nonlinear 
simulation and linear time-varying covariance propagation programs. 


* 
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Table 1. SPIFR Experimental Matrix 
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VALIDATION Of UH-1H APPROACH AND 
LANDING SIMULATION PROGRAMS 


A first step In the pursuit of an analytical methodology for de- 
termining pilot acceptance Is the development and validation of offline 
computer simulation programs capable of mimicking automatic as well as 
manual flight director approach and landing controls strategies and 
closed-loop pilot-aircraft system performance. The use of two existing 
UH-1H approach and landing simulation programs resident on the NASA- 
Ames IBM 360 computer was stipulated at the Initiation of this research 
effort. They are: 

1) A UH-1H nonlinear simulation program for the Monte-Carlo 
simulation of the 'pproach and landing task - hereafter referred 
to as the Nonlinear Covariance Computation (NCC) program; and 

2) a linear time varying simulation program for the propagation 
of approach path dispersions (l.e., mean and standard devia- 
tion) during the approach and landing task - hereafter called 
the Linear Time-Varying Covariance Progagatlon (LTVCP) program. 

Therefore, the first two tasks In this study effort were to validate 
the approach performance prediction capability of the existing NCC 
and LTVCP simulation programs using data from the previously completed 
SPIFR experiments. The NCC program validation vj discussed first followed 
by an analysis of the LTVCP program. 


Nonlinear Covariance Computation (NCC) 

Program Validation 

Description . The NCC program provides a closed loop simulation of the 
UH-1H approach and landing task. The 6 degree of freedom aircraft equa- 
tions of motion for the UH-1H as given In reference 7 are used to simulate 
the basic helicopter nonlinear dynamics and are incorporated into a larger 
and general IBM 360 program existing at NASA/ Ames, called "UNCLE". Four 


different control options are provided corresponding to*. 

1) Basic helicopter (WSB) with 3-cue flight director and raw 
data for the directional axis (BH/FD), 

2) Basic helicopter (WSB) with CSS, 3 cue flight director and 
automatic directional axis (CSS/FO), 

3) Basic helicopter with CSS, 2 cue flight dlrjcor, and automatic 
vertical and direction 1 axis (CSS/FO/AC), 

4) Basic helicopter w 4 th CSS, 3 cue flight director and fully 
automated system (CSS/Auto). 

The control stick steering configuration consists of pitch, roll and ver- 
tical SCAS and an automated directional control system. Block diagrams 
for the longitudinal, vertical, lateral and directional systems are 
shown In Figures 7-10, respectively. Details on the design of the CSS 
augmentation and flight director may be found In reference 3. 

A velocity command, V c , corresponding to the constant attitude de- 
celeration profile Is simulated to allow for decision heights of 100 ft, 

50 ft and 0 ft respectively. The velocity command is generated by solving 
an Implicit algebraic equation between command velocity and distance to 
go. 


A Dryden mode, lias shown In Figure 11 is used to simulate the three 
orthogonal turbulence components Ug, Vg and Wg In aircraft body axes. A 
root mean square value of 1.7 ft/s Is used for each of the three wind gust 
standard deviations. Mean winds or shears were not implemented in the off- 
line program. 

A plotting routine (PLOTTA/DISSPL) is used to allow plotting SPIFR 
data from digital tape and NCC program outputs on the same page (e.g., 

Beta or UCC plots). This permits direct visual evaluation of the quality 
of fit between NCC program outputs and piloted simulation data. 
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Fiyure 8. Vertical System 
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Figure 9b. Automatic and Flight Director Guidance 
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Figure 10. Directional System 
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Figure 11. Dryden Turbulence Models 


Some of the above options were not available on the existing IBM 360 
program and were added to make realistic comparison to experimental data 
feasible. The NCC program was then exercised to match experimental data 
in the automatic mode and the CSS/FD mode. A slngle-input/slngle-output 
describing function model of the form shown in Figure 12 was implemented 
in each of the three flight director axes. In this model, k is the pilot 
gain, -z is the zero, -p is the pole and x is the effective time delay. 
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i = pitch, roll or vertical axis 
Figure 12. Input/Output Describing Function Pilot Model 


The selection of this simple model was merely for the purpose of 
initial verification of the flight director mode. More realistic re- 
presentation of the pilot's decision making and control behavior would 
be essential in order to reproduce some of the finer details of human 
operator response data. 
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However, the existing NCC program has several problems associated with 

It; the most significant being the lack of adequate documentation for the 

UNCLE program that Is used to represent tho UH-1H nonlinear equations of 

£ 

motion. This made It very difficult to comprehend, analyze or modify over- 
all NCC program operations. In addition, the program Is Inefficient for 
use In the Monte Carlo mode. For example, a single 120 s run using the NCC 
program on the IBM 360 takes, approximately, 170 CPU seconds (CSS/Auto mode) 
to complete at an Integration step size DTP of 0.1 s. The flight director 
mode takes 10 ~ 20 per cent longer. An ensemble calculation using 10 runs 
would take approximately 30 minutes CPU time on the IBM 360. This amount 
of time Is prohibitive for routine use. Two options were considered: (a) 
convert the existing UH-1H nonlinear simulation program (UNCLE version of 
the NCC) to the CDC 7600, or (b) start from scratch and develop a separate 
IBM 360 UH-1H nonlinear simulation program using only those parts of the 
UNCLE program that are useful for this study. The first option was not 
seriously entertained in view of the dearth of documentation on the existing 
UNCLE program and the inefficient nature of the overall code. Instead, it 
was decided to develop a separate NCC program on the IBM 360 computer using 
only those elements of the overall UNCLE program (l.e., the UH-1H equations 
of motion) that were needed for the intended simulation. This program was 
developed and response time history results are indistinguishable from the 
original NCC outputs. However, the new version takes 35 CPU seconds at a 
step size DTP or 0.1 s, for a single 120 s CSS/Auto run. This is at least 
five times more efficient than the original version of the NCC program for 
the IBM 360 computer. Additional cost savings may be obtained by conver- 
sion to the CDC 7600 computer. 

Simulation Results 

The new version of the NCC program was exercised for the CSS/Auto, 
CSS/FD and CSS/FD/AC configurations. Figures 13 - 14 show a comparison 
of the NCC program outputs ( i . e . , ensemble mean and standard deviation) 
with experimental data for the CSS/Auto/DH = 0 ft mode for the no wind 
(IWIND = 0) and wind ( I WIND = 1) conditions. A total of 11 flight 
variables are displayed with standard nomenclature used to label each 
output. The following observations are worth noting: 

1. The mean response comparisons show an excellent fit between 
the NCC program outputs and the SPIFR data. 
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Figure 13. NCC Program Simulation Results: CSS/AP/NO WIND/DH = 0 ft, 10 runs 
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Figure 13. NCC Program Simulation Results: CSS/AP/NO WIND/DH = 0 ft, 10 runs 
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Figure 14. NCC Program Simulation Results: NCC/AP/WIND/DH = 0 ft, 10 runs 
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Figure 14. NCC Program Simulation Results: CSS/AP/WIND/DH - 0 ft, 10 runs 
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Figure 14. NCC Program Simulation Results: CSS/AP/WIND/DH * 0 ft, 10 runs 
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2. The reference heading angle ij/ ref . for the SPIFR simulation 
experiments was >7.3 degrees (corresponding to Crow's landing 
runway); hence the apparent offset error In the PSI responses. 

3. The NCC program does not Include the blending logic for cap- 
ture and flare used In the SPIFR experiments. This discre- 
pancy Is responsible for the differences between the two out- 
puts at glide slope capture and Initiation of deceleration. 
However, this omission has no significant Impact on the overall 
time response history. 

4. The automatic mode ensemble standard deviation data from SPIFR 
experiments Is significant even In the absence of wind gust 
disturbances as evidenced In Figure 13 c,d. In contrast the 
NCC program output standard deviations are zero as would be 
expected. The non-zero standard deviations must be the result 
of the peculiar hybrid mechanization of the simulation In the 
SPIFR experiments. Run to run variability and hence the non- 
zero standard deviations may be due to (1) the lack of synchro- 
nization between the on-board 1819B Sperry computer and the 
EAI 8400 machine used for simulating the aircraft dynamics; and 
(2) the variability in the Initial conditions at the start of 
each run. Other unlikely sources irclude the quantization and 
truncation errors in the on-board computer or a variable time 
delay In the A/D and 0/A mechanizations. 

5. Comparison of the aul xnatic mode standard deviation results 
with and without wind (Figures 14 c,d and 13 c,d, respectively) 
shows that the relative contribution of wind disturbances Is 
roughly equivalent to that caused by the specific V/STOLAND 
system digital computer implementation. 

6. Standard deviations of the ground speed , VG, and pitch attitude, 
THETA show a sudden decrease after initiation of flare at ap- 
proximately 4000 ft to go. This is consistent with the inple- 
mentation of the velocity guidance law where the velocity feed- 
back used for comparison against commanded speed V c switches 
from airspeed VRW to ground speed VG. Hence groundspeed VG 
shows a higher standard deviation before the onset of decelera- 
tion. 
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Results for three CSS/FO approaches are shown In Figures 15-17. The 
zero decision height response comparison Is shown In Figure 15, where the 
experimental data Is for pilot number 1, the only pilot who attempted the 
zero/zero approach to a hover. Comparisons for the decision height of 50 ft. 
are shown In Figure 16. The 2 cue flight director results corresponding 
to the automatic collective configuration are presented In Figure 17. SPIFR 
data for pilot number 2 Is displayed In both Figures 16-17. The following 
characteristics may be noted: 

1. The mean response comparisons In the flight director mode are 
satisfactory In all axes. The zero/zero case shown In Figure 15 
shows the best match to the experimental data. The apparent 
error In the PSI response fit Is due to the fact that the re- 
ference heading angle for the SPIFR simulation was -7.3 
degrees (Crows Landing runway) while It Is set to 0 degrees 

In the NCC program. 

2. The pitch axis mean response comparisons show significant 
differences during the deceleration phase. The deceleration 
command Is provided tr the pilot through the pitch flight 
director bar PFDIN. At the Initiation of flare (as Indicated 
by a green FLR light on the ADI), the pilot must respond by 
attempting to track the pitch flight director signal while 
keeping the aircraft on the prescribed flight path. Any 
delay in responding to the pitch flight director command re- 
sults In larger command attitude to compensate for the loss 
of distance remaining to decelerate to a hover. Pitch axis 
data In Figures 15-17 demonstrate this type of behavior on 
behalf of the human pilot, e.g., see THETA, THCOM, PFDIN and VG. 

A pure gain (no lead/lag or time delay) model for the pilot 

Is assumed in all three control axes. Improvements in fit may 
be obtained using a more realistic pilot model. 

3. Response differences in the lateral axis variables. In par- 
ticular, the localizer error DY and bank angle PHI, are due to 
the high variability of the initial conditions evident in the 
SPIFR simulation data. In the SPIFR experiments, the automatic 
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Figure 15. NCC Program Simulation Results: CSS/FD/WIND/DH = 0 ft, 10 runs 
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Figure 16. NCC Program Simulation Results: CSS/FD/WIND/DH = 50 ft, 12 runs 


44 





NCC Program Output 
SPIFR Data; Pilot No. 2 




Figure 17. NCC Program Simulation Results: CSS/FD/AC/WIND/DH * 50 ft, 12 runs 
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Figure 17. NCC Program Simulation Results: CSS/FD/AC/WIND/DH * 50 ft, 12 runs 
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mode was used for the first few seconds In the run before the 
pilot took over manual control. Therefore the variability of 
the Initial conditions Is partially the result of the run-to- 
run variations In engaging the CSS/FD mode. In addition, data 
Is recorded from a distance of 10,200 ft from the landing pad 
while the actual runs were Initiated at a distance of 13,200 ft. 

As a result Initial conditions on DY when data recording was 
initiated are typical large (non zero) and display a large 
standard deviation as witnessed by the data In Figures 15-17. 

NCC program initial conditions for each run In the Monte Carlo 
process are set to zero or equivalent trim conditions. A fair 
comparison demands, wherever feasible. Identical factors In the 
two simulations. If this were possible, errors between the two 
sets of data would be due to the discrepancies In the modeling 
of the pilot control strategy alone, and efforts could be 
directed towards developing a more realistic representation. 

4. The standard deviation results display significant differences 
between the SPIFR data and the NCC program outputs, the latter 
being highly optimistic. As mentioned earlier in discussing 
the automatic mode results, some of the differences are due to 
the V/STOLAND digital implementation characteristics which are 
not modeled in the NCC program. However, a major source for 
discrepancy lies in the highly idealized representation of the 
human pilot as a pure gain servomechanism. Clearly, the human 
pilot does not behave as a robot or an autopilot. A more 
realistic model for the pilot would include an adequate re- 
presentation of human information acquisition (scanning strategy), 
data processing and decision making capabilities. Such a model 
would show signifigant run-to-run variability and more closely 
match SPIFR experimental data. 

Results presented here demonstrate the validity of the NCC program. 
Observed discrepancies have a rational explanation and can be alleviated 
by improving some of the Individual modules in the overall program. Within 
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the scope of this effort, an Improved representation of pilot response be- 
havior Is considered more Important than developing more realistic models 
for the V/STOLAND system Implementation, Consequently, the NCC program 
Improvements are directed towards development of a realistic pilot model 
for use In formulating an analytical methodology for determining helicopter 
precision approach requirements. 


Linear Time-Varying Covariance Propagation 
(LTVCP) Program 

An investigation of the existing LTVCP IBM/360 program was con- 
ducted with the goal of determining its usefulness for application to the 
SPIFR simulation data. The existing LTVCP program is suitable for solving 
the standard covariance propagation equations as follows: 

Given the linearized aircraft equations of motion 

x * Fx + Gu + Tw (8) 

where x = perturbed state vector, u = input, w = N(p,Q), the program i ves 
the equations for the mean x - E(x) 

t 

x = FiT + Gu + Tp (9) 

and covariance X = E(x-x) (x-7) T 

X = fx + xf t + rqr (10) 

In general F,G, and r are matrices obtained by linearizing the aircraft 
nonlinear equations of motion around some operating point (V Q , y) where 
V. is airspeed and y is the flight path angle. Thus, 

u 

f, G, r = f(v a , y) (id 

The existing LTVCP program is not suitable for application to the SPIFR 
experimental scenario because: 

1) The program only selves for the frozen point or steady state 
solution for the state covariance using 
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x ■ o ■ fx + xf t ♦ rqr T 


(12) 


where F and r are obtained by Interpolation between the nominal 
flight conditions - namely trim values of airspeed V fl , flight 
path angle y and longitudinal deceleration Vj. Steady state or 
frozen point covariance computation Is acceptable If changes In 
nominal operating flight conditions are slow compared to the 
system response time. For the SPIFR simulation, the changes In 
the commanded nominal flight conditions at glideslope capture 
(y « 0 to y * -6*) and Initiation of the constant attitude de- 
celeration maneuver (0 nom ■ 3.94° to 6 nom s T*) occur too 
rapidly compared to the system response time to complete the 
transient maneuvers. For example, approximately 6 seconds are 
required to maneuver the pitch attitude of the aircraft from 
the trim value of approximately 3.94° during the glide slope 
trajectory phase to the conmanded value of 7° following the on- 
set of the CADP phase. Therefore, frozen point solutions are 
not applicable to the analysis of the precision approach and 
landing task. 

2. The program allows for interpolation of the system matrices 
in only one independent variable - only V. in this applica- 

Q 

tion. This is unacceptable for the SPIFR simulation reference 
flight trajectory where changes in flight path angle and 
airspeed may occur, independently or simultaneously. 

3. The program has a peculiar structure which requires a specific 
input-output format. This puts undue restrictions on the prac- 
tical use of the program for application to the SPIFR simula- 
tion structure. 

4. The program by definition requires the use of data where time 
is the independent variable. Time ensemble data must, there- 
fore, be used for validation. However, ensemble averaging with 
time as an independent variable is not appropriate in situa- 
tions where the guidance commands are computed as a function of 
distance to go. In the SPIFR simulation, velocity and attitude 
commands are based on distance to go from the hover point. This 
could lead to situations where successive runs may be identical 
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with respect to distance to go as the Independent variable but 
could be advanced or delayed with respect to each other In time. 
Time ensemble statistics would clearly be Inappropriate In such 
circumstances. 

5. The program only provides the ensemble mean x p and covariance x p 
of the perturbed state vector x p . The total system response x 
Is the sum of the nominal x n and perturbed state x p . Thus, 


x(t) 

■ *„(t) ♦ Xptt) 


and x(t) 

■ x „(‘> ♦ yo 

(13) 

X(t) 

* y 1 * 



Thus knowledge of x n (t), the nominal state vector is essential 
for computing the total system vector. However, x n (t) can only 
be obtained by solution of the complete nonlinear simulation of 
the system as In the NCC program. The principal justification 
of using linear time varying covariance propagation is undermined 
if a nonlinear simulation is required for computing the nominal 
response. 

6. Finally, the program also gives highly optimistic results for 
the state covariance propagation. Actual state vector disper- 
sions are caused by wind as well as other factors such as the 
quantization and truncation errors due to the particular on- 
board computer implementation. The available linearized models 
do not allow for error models for digital computer error sources. 

In conclusion, comput ion of approach and landing performance using 
linear time varying covariance propagation methods is not recommended. The 
problems with the approach are fundamental and not program specific. Specific 
limitations (1) - (4) given above can be remedied by developing a new pro- 
gram which uses distance as the independent variable for covariance propaga- 
tion. However, restrictions (5) and (6) above are more fundamental to the 
methodology itself and cannot be removed creating a new program. 
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Linear time varying covariance propagation programs apply as long as 
the perturbations In the system state variable stay reasonably close to the 
nominal trim values. Application to the decelerating approach to a hover 
task Is complicated by the requirement of large sudden changes In the opera 
ting point or helicopter trim conditions during glide slope capture and 
Initiation of flare. The resulting dynamic trim change cannot be described 
as a string of frozen point trim (static or dynamic) data, and must be 
obtained through a nonlinear system simulation. Therefore, an approach 
using a nonlinear simulation of the complete system as In the NCC program 
Is desirable, especially for systems Involving dynamic maneuvers such as 
reference trajectory capture and deceleration to a hover. 
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PILOT ACCEPTANCE ANALYTICAL METHODOLOGY 


The overall approach to the formulation of an analytical methodolgy 
for determining pilot acceptance consists of the following four steps: 

1. Develop and validate an off-line simulation of the helicopter 
approach to a hover task In the automatic and flight director 
modes, 

2. Select a generic pilot model structure that Is applicable to 
manual flight director or raw data modes of operation, 

3. Formulate candidate pilot model based workload metrics and 
evaluate their adequacy through simulation of the performance- 
workload tradeoffs, and 

4. Summarize the results of the study by formulating a plausible 
analytical procedure for determining the pilot acceptance 
ratings of competing helicopter guidance, control and display 
systems. 

Based on the results of the last section, the NCC program has been 
sufficiently validated to serve as the benchmark simulation of the heli- 
copter precision approach to a hover task. 

The steps (2)-(4) are based on the implicit assumption that pilot 
acceptance is related to the information, decision and control require- 
ments imposed upon the human pilot in satisfactorily accomplishing a given 
task. In other words, determination of pilot acceptance or opinion rating 
is assumed to depend upon two factors: 1) a measure of system performance 
as perceived by the human pilot, and 2) a metric for human information, 
decision and control requirements usually lumped into the undefined but 
ubiquitous word "workload". However, both system performance and pilot 
workload depend upon the open loop aircraft characteristics and the struc- 
ture of the information processing, decision making and control strategy 
adopted by the human pilot. Thus, the need for a generic dynamical model 
for the pilot is apparent. Assuming this Is feasible, pilot opinion 
rating may be specified as a linear or nonlinear function of the system 
performance measures and pilot model parameters. 
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The pilot model -based approach has been proposed by a number of 
Investigators using a variety of pilot model structures. For example, 
the so-called "Paper Pilot" (9] approach, Is based upon the use of des- 
cribing function pilot model and assumes that the pilot opinion rating 
R Is related to performance J and workload W by a functional of the form 

R - J + W + 1 

where system performance J find workload W represent multlattrlbute utility 
functions which depend on the Individual components of a performance 
vector s and pilot model parameter vector e. Figure 18 shows an application 
of the paper pilot concept to a longitudinal precision hover task [10]. 

In this application, system performance J is decomposed as a function of 
the rms values of longitudinal deviation (o x ) and pitch rate (c^), where 
as pilot workload W is defined to be a nonlinear function of the pilot's 
pitch attitude lead (T^ 0 ) and longitudinal position error lead (T LX ) para- 
meters. 

Pilot workload metrics and pilot-rating functionals in terms of the 
optimal control pilot model have also been proposed [11,12]. For example, 
in one approach proposed by Levi son [12] pilot workload is assumed to be 
proportional to the fraction of attention devoted to the flying task as a 
whole. A visual curve fit or regression analysis is used to find the par- 
ticular form of the functional relating pilot opinion rating to the system 
performance and workload metrics. 

Another method based on the optimal control pilot model [13] suggests 
a pilot opinion rating functional that is monotonically related to the 
quadratic index of performance used In the model formulation. Thus, 

R = U( I ) 

t T 

where I = e!}™| J (vV + u T Ru)dt 
' o 

where y (t ) represents the pertinent linear combination of the aircraft 
states and u represents the pilot's control vector (aircraft control 


(14) 

(15) 
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Inputs). Note, that Implicit In the definition of the performance Index 
I Is the partition Into the two components of system performance J and 
pilot workload W as 


0 ■ Efi f (y T 0y)dtJ 

0 

and 

W ■ e [}!2 y T ( u T Ru)dtJ 

0 


( 16 ) 


( 17 ) 


These model based pilot rating functionals are usually derived for 
stationary or frozen point conditions using a limited set of data. They 
merely reflect the existence of a correlation between actual pilot opinion 
ratings and model -based values, and should not be construed to have pre- 
diction capability. The predictive capability of these model -based tech- 
niques is intimately dependent upon the validity of the pilot model struc- 
ture used, and the assumed values of the pertinent parameters. 


The true value of the model based approach towards determining pilot 
acceptance lies in the formulation of an acceptance methodology. The 
weakest link is the process of selecting an appropriate pilot model for 
the task under evaluation. Helicopter approach to a hover under IMC 
is a very complex and demanding piloting task requiring deterministic 
command maneuvers such as localizer and glideslope capture, and constant 
attitude deceleration along a prescribed flight path. 


Existing pilot models - namely the quasilinear describing function and 
the standard optimal control models are not suitable for application to 
the helicopter approach to hover task. These models have been developed 
primarily for description of human input-output behavioral characteristics 
in an aggregate or ensemble sense and are not suited for modeling response 
to deterministic sample Inputs on an individual run basis. Therefore, the 
following structure for a generic dynamical model for the pilot that in- 
corporates decision logic for attention allocation and control involvement 
is formulated and incorporated as a module into the nonlinear covariance 
computation (NCC) program. Description of the proposed model structure is 
given next. 
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Pilot Model Formulation 


The human pilot In a realistic flight control environment must con- 
tlnually perform the following three tasks: 1) lnformatlng acquisition, 

2) Information processing, and 3) aircraft control. Information acquisi- 
tion Involves the visual scanning of the cockpit displays to collect raw 
data necessary for the particular flight segment. This raw data are then 
processed by the human brain to generate an estimate of the significant 
aircraft states. Finally, the estimated states may require further 
aggregation In the form of a decision function, or decision rule, upon 
which decisions regarding selection and Implementation of a control 
strategy are based. A pilot model Incorporating these features can be 
formulated as shown In the functional block diagram representation of 
Figure 19. The selection of a specific form or structure for each of 
the four block diagram components Is assumed to be made at a highet 
cognitive level (the brain) called the decision making center in Figure 
19. 


The pilot model formulation of Figure 19 is an attempt towards 
homomorphic representation of human information acquisition, processing and 
decision making behavior. This approach is significantly different from 
currently used pilot models based on describing function methods and 
optimal control theory. However, it is felt, that a functionally homo- 
morphic model as formulated in Figure 19 is essential for analysis of 
dynamic, nonlinear time varying tasks such as the helicopter decelerating 
approach to a hover. Similar modeling concepts have been proposed 
by investigators In Europe [14,15] and the United States [16]. The 
function of the individual model components is described in the following. 


Scanning Rule The human pilot acquires primary information about the 
state of the aircraft from the cockpit displays, through visual scanning 
of the individual Instruments. Several models for the human scanning 
strategy have been proposed in the literature [17-20]. The models are 
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Figure 19. Functional Model for the Human Pilot 


based on information- theoretic concepts and represent various alternative 
hypotheses for human visual sampling behavior. However, formulation of a 
scanning rule Is an Intrinsic part of human pilot behavior, and, therefore, 
must be included in any realistic pilot model. 

State Estimator The raw data y m obtained through scanning the aircraft 
instruments must be processed to estimate a state vector z which may con- 
sist of display variables as well as additional states such as the air- 
craft angular velocities not directly displayed on the cockpit panel. 

Decision Function The decision function DF consists of an aggregation of 
the estimated state variables in z, and is used in the selection of a con- 
trol strategy for the given task. The human pilot has a choice of applying 
control in one or more axes, either simultaneously or sequentially. The 
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decision function reflects the degree of urgency or seriousness of each 
control axis In terms of Its r*eed for pilot corrective action or control 
Inputs. Decision functions have otherwise been referred to as urgency 
functions by Onstott [16] or Instantaneous seriousness functions by Cavalll 
[15]. 

C ontrol Strategy Several control strategies are available tc the human 
pilot, ranging from slmulataneous , continuous or Intermittent, Involvement 
In all control axes, to sequential and Infrequent control based on the 
urgency for corrective action. 

Decision Making Center The purpose of the decision center Is to perform 
the high level supervisory task of coordinating the Individual component 
functions In an efficient manner. The decision center has the authority 
to modify the scanning rule, the state estimator structure, the decision 
function formulation and the control strategy selection criteria depending 
upon the overall task requirements. 

A framework for describing human pilot Information processing, con- 
trol and decision making response characteristics has been proposed. Utility 
of this approach to the determination of system performance, pilot work- 
load and eventually pilot acceptance of a helicopter precision approach 
to a hover task Is considered In the next section. 


Application to the Helicopter 
Approach and Landing Task 


The purpose of this section is to demonstrate the applicability of 
the generic pilot model structure proposed above to a realistic flight con- 
trol task as evidenced in a helicopter decelerated approach to a hover under 
IMC. The CSS/FD mode In the SPIFR simulation experiments Is used to Illus- 
trate the key concepts. 
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The helicopter pilot uses Information from the available cockpit dis- 
plays to generate the appropriate control conmands to the aircrafts so as 
to accomplish the desired flight control task objectives. In the CSS/FO 
mode, the pilot has three flight director Inputs corresponding to the pitch, 
roll and collective axes. For the purpose of this study, the pilot model 
Is restricted to slngle-lnput/slngle-output representations for each of the 
three control axes. The general structure of the single-axis pilot models 
Is shown In Figure 20. where (Xj.Uj) 1 ■ 1,?,3 represent the three equiva- 
lent flight director Inputs and corresponding pilot command outputs for the 
pitch, roll and collective axes, respectively. 

In the raw data mode, the three equivalent Inputs x^ must be treated 
as some linear or non-linear combination of the appropriate raw data signals. 
For example the roll and vertical axes equivalent Inputs Xg and x 3 In the raw 
data mode could easily be justified as the raw localizer and glide slope 
errors or filtered through some form of a dynamic compensator (e.g., lead 
network). The pitch axis equivalent Input Xj Is more difficult to hypoth- 
slze. However, It may be reasonable to assume that the human pilot has de- 
veloped some Internal model for the desired velocity profile as a function 
of the radio altitude. This may be justified because the pilot does not 
have ground speed (or range rate) or range from hover pad available In the 
SPIFR simulation and must therefore rely on some form of heuristic relation- 
ship between airspeed and radio altitude h (above hover pad). Errors 
must be expected In the presence of wind and due to non-zero glide slope 
tracking errors. An equivalent fourth directional flight director signal 
must be synthesized for the basic helicopter mods. 

The model allows for two Independent periodic samplers and S c - one 
on the display variables called the display sampler (S d ) and another on the 
pilots control output (S ). The human pilot does not necessarily apply 
corrective control Inputs to the aircraft at the same rate at which Infor- 
mation from the display variables Is being obtained. Thus, control sampling 
rates reflecting degree of pilot involvement In the control task need not 
be Identical, and may inHfact be lower than display sampling rates. For 
mathematical convenience, the two samplers S c and S d are as'umed to operate 
at sampling Intervals T c and T d , respectively, such that T c = (Integer)^. 

The subsequent analyses are based on the assumption that pilot workload is 
related to the degree of pilot involvement in the information acquisition, 
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Figure 20. Single Axis Pilot Model Structure 
processing and control tasks. 

Dynamic compensators In the form of the classical describing function 
model (gain, lead/lag and time delay) are allowed between the pilot's com- 
mand outputs u^ and aircraft control Inputs. 

Several options are available to the supervisory decision maker In 
selecting the appropriate control strategy as described below. 


Strategy Option 0: (ISTRAT * 0) Independently adjustable and 

periodic control sampling rates 

This option permits simultaneous control of all three axes with 
different levels of control involvement. Thus, if DTP A Integration step 
size, then the level of control Involvement as defined by the control sampling 
period T c ( 1 . e . , sampling Intervals at which display variable is used for 
control purposes) in the pitch roll ana vertical axes may be selected to 
be (INPTCH*DTP), ( INR0LL*DTP) and INC0LL*DTP) seconds, respectively, 
where the three multipliers INPTCH, INROLL and INCOLL are constrained to 
be positive integers. In this option, the estimator and decision func- 
tion element are absent. 

The maximum control sampling rate, f = i • for each axis Is (1/DPT) Hz. 

c i c 

Typically, the pilot aircraft system tends to become less stable with de- 
creasing sampling frequency. Some adjustment, usually a reduction in the 
corresponding pilot gain K is required to increase the stability boundaries. 
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Strategy Option 1: (ISTRAT ■ 1) Single channel hypothesis with 

' state dependent control sampling rates 


This strategy Is based upon the assumption that the human pilot can 
only operate one control at any given time with the others In hold or 
slowly decaying towards trim condition. A higher level decision logic 
Is assumed for switching between the 3 axes and a possible fourth mode 
corresponding to an auxiliary task or any other activity. The decision 
logic assumed decision functions DF^ (i * 0, 1,2,3) for each axis as 
follows (from Onstott - reference 16): 

DF i = |a 1^1 + B sgn(x 1 )x i | , 1 « 1,2,3 (18) 

* UMAX , 1=0 (19) 

The allocation rule Is to choose axis i for active control if 

DF. = Max |DF . | (20) 

1 j = 0,1, 2, 3 3 

with the remaining control outputs in hold or exponentially decaying to 
trim with some time constant. Only the collective controls is put in 
hold; the pitch and roll cyclic control are set to decay towards their 
trim values. The decision function DFq = UMAX may be used as a measure 
for the reserve attention capacity of the pilot. If UMAX = 0, attention 
is shared, one axis at a time, between the pitch, roll and vectical con- 
trol axes, depending upon which o* the 3 decision functions DF.. is largest. 

A very large value for UMAX would result in essentially open loop control 
of the aircraft and would lead to instability. The largest value of UMAX 
allowable within allowable performance constraints on all 3 axis serves 
as a measure of reserve attention capacity and hence pilot workload. 

The decision function computation and comparison may be carried out 
every sample (DTP seconds) or every (INUMB2*DTP) seconds where INUMB2 can 
be set to any positive intecer value. 
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An option Is provided whereby the continuous control Involvement In 
any one axis Is limited to a prescribed time or consecutive number of 
samples INUMB3 after which the control must be switched to one of the re- 
maining alternatives according to the same decision function rule. 

In addition control sampling rates reflecting degree of pilot In- 
volvement In the axis under active control may also be adjusted by mani- 
pulating the appropriate multipliers INPTCH, INROLL and INCOLL. 


Strategy Option 2 : (ISTRAT * 2) Simultaneous control of all axis 

with state dependent control sampling rates 

This control strategy Is akin to control strategy 1 with the exception 
that all three control axes may operate Independently. The decision func- 
tions DF. are defined as in option 1. The attention avocation rule 
operates on each ixis independently of the others. Decision to control 
axis i is made if 

DF. > ULIM. i = 1,2,3 (21) 

Thus, at any given time none, one, two or all three axes may be under 
active control. The limits ULIM.. represent threshold values for deciding 
between active control of axis i or open loop control. 

As in option 1, the periodic interval at which the decision functions 
are calculated and compared to their respective thresholds is adjustable 
by varying the positive integer INUMB2. 


Candidate Pilot Workload Metrics 

Pilot workload may be assumed to depend upon the minimum required level 
of pilot involvement in the information acquisition and control actuation 
tasks fos a given acceptable level of system performance. Three different 
hypotheses or options have been presented as candidates for describing 
human pilot input-output characteristics. Consequently, three different 
workload metric vectors may be defined corresponding to the three options 
as follows. 
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Option 0: Here the model parameters INPTCH, INROLL and INCOLL specify 
the control sampling intervals In the pitch, roll and vertical axis. 

Pilot workload In the pitch axis may be assumed to be Inversely re- 
lated to the maxi mum allowable sampling Interval INPTCH for a given 
level of longitudinal system performance. Similar relationships apply 
to the roll and vertical axis. Thus workload metrics for each c. 
the three control axes may be defined as ' 

WLPTCH - hj (INPTCH) 

WLROLL » h 2 (INROLL) (22) 

WLCOLL = h 3 (INCOLL) 

where hj, h 2 and h 3 are Inverst monotonic functions of their respec- 
tive arguments. Total pilot workload WL may be defined as some linear 
or nonlinear function of the individual workload metrics; thus 

WL = h(h 1 ,h 2 ,h 3 ) (23) 

Option 1: This option is characterized by decision functions DF 4 , 

i = 1,2,3 and the reserve attention capacity measure DF r = UMAX. 

Thus, UMAX may be used as a workload index for the overall task. Re- 
lative workload of the three control axes may be defined as in Option 0 
in terms of the relative fractions of the total time spent on each 
axis. These numbers must be obtained through Monte Carlo simulation 
of the NCC program. 

Option 2: This option allows for simultaneous control of one or more 

axes depending upon the values chosen for the three decision function 
limits ULIM.., i = 1,2,3. Pilot workload may again be defined in terms 
of the relative attention requirements of each axis for a given level 
of system performance. These fractions of total time spent on each 
axis must be obtained through Monte Carlo computer simulation of the 
NCC program. 
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Model Simulation Results 


The generic dynamical model for the pilot as described above was In- 
corporated as a module Into the nonlinear covariance computation (NCC) pro 
gram. A program description In the form of a flowchart Is given In the 
Appendix. A flowchart and listing of subroutine F2T Incorporating the 
pilot decision/control logic subroutine DECSN Is also Included. The three 
different decision and control strategy options corresponding to ISTRAT 
of 0, 1 and 2 were simulated In the CSS/fllght director configuration to 
study the sensitivity of the system performance to key control strategy 
parameters. The different cases Investigated for each mode are described 
bel ow : 

1. Strategy Option ISTRAT = 0 

The following sets of parameters were simulated. 


Set No. 

INPTCH 

INROLL 

INCOLL 

1 

1 

1 

1 

2 

5 

5 

5 

3 

10 

10 

10 

4 

1 

5 

10 

5 

1 

10 

10 

6 

5 

10 

10 


Figures 21-26 show the ensemble mean and standard deviation responses for 
11 key system variables. 

The results demonstrate the sensitivity of the system variables to 
the attention allocation parameters INPTCH, INROLL and INCOLL for the pitch, 
roll and vertical axes, respectively. The ensemble plots are terminated 
around the time that the system goes unstable. The parameters INPTCh 
INROLL and INCOLL provide a measure of the degree of control involvement 
in the pitch, roll and vertical axes, respectively. Pilot workload may 
be assumed to be related to the degree of control involvement required 
to achieve a prescribed level of performance. Decreasing the control 
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Figure 21. NCC Program Simulation Results: CSS/FD/WIND/DH = 0 ft 
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Figure 22. NCC Program Simulation Results: CSS/FD/WIND/DH = 0 ft 
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Figure 73. NCC Program Simulation Results: CSS/FO/WIND/DH * 0 ft 
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Figure 23. NCC Program Simulation Results: CSS/FD/WIND/DH = 0 ft 
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Figure 23 . NCC Program Simulation Results: CSS/FD/WIND/DH 
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Figure 24. NCC Program Simulation Results: CSS/FD/WIND/DH s 0 ft 
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Figure 25. NCC Program Simulation Results: CSS/FD/WIND/DH 
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Figure 25. NCC Program Simulation Results: CSS/FD/WINU/DH = 0 ft 
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Figure 26. NCC Program Simulation Results: CSS/FD/WIND/DH = 0 ft 
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Figure 26. NCC Program Simulation Results: CSS/FD/WIND/DH * 0 ft 
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Involvement equally In each a;t1s by simultaneously Increasing INPTCH, 

INROLL and INCOLL from 1 through 10 results In early Instability as seen 
1r Figures 21-24. The relative difficulty of the Individual control 
axes may be measured by comparison of the maximum allowable values of 
INPTCH, INROLL and INCOLL compatible with acceptable overall system per- 
formance. Results show that the pitch axis Is the most sensitive to the 
degree of attention allocation and control Involvement as evidenced by 
the rapid degradation of system performance with Increasing values of 
INPTCH. 

Additional single run (ICARLO * 1) simulations for varying values of 
the three parameters show that the limiting values of INPTCH, INROLL and 
INCOLL are 7, 15 and 20 respectively. 

Eventually the values of MPTCH, INROLL and INCOLL In addition to 
being Independent of each other may be allowed to vary with the nature 
of the control task or as a function of the distance to go from the hover 
pad. Thus, the 3 parameters could vary from relatively large values during 
the cruise phase towards smaller values (approaching 1 } as the flight 
progresses through glldeslope capture, deceleration and flare to hover 
over the prescribed landing pad. This would be compatible with the In- 
creasing difficulty or workload level of the piloting task with decreasing 
range from the hover pad. Eventually, the 3 parameters may be programmed 
to be adaptive to system performance according to some adaptive control 
law mimicking human adaptive response capabilities. 

2. Strategy Option ISTRAT ■ 2 

In this option, the degree of control Involvement In each axis Is 
assumed to be directly related to the magnitude of the Instantaneous errors 
In following the pitch, roll and vertical flight director commands ( i ,e . , 

Oj * ag * 013 * 1, 0 j * $ 2 s £3 “ 0 in Eq. 18). The following parameter 
sets were simulated to illustrate the effects of varying the performance 
threshold levels ULIM(I), I * 5,2,3. 


£- 2 . 
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Set No. 

ULIM(l) 

UL!M(2) 

ULIM(3) 


(l"-15*j 

(1"»66 # ) 

(l-20ft/s) 

1 

0.0 In 

0.0 In 

0.0 In 

2 

0.034 In 

0.03 In 

0.05 In 

3 

0.067 In 

II 

ll 

4 

0.134 In 

II 

II 


Figures 27-29 show the simulation results for set numbers 2-4 above. Set 
1 results are Identical to those In Figure 1 for ISTRAT ■ 0, INPTCH ■ INROLL 
s INCOLL * 1 as should be expected. Larger threshold values result In 
Increasingly oscillatory behavior In the corresponding control axes. This 
Is evidenced In Figures 27-29 which demonstrate Increasing oscillatory be- 
havior In the pitch axis as the corresponding pitch flight director thres- 
hold parameter ULIM(l) Is Increased from 0.0 In (equivalent to continuous 
control) to 0.134 In (or a deadzone of +2° In responding to the pitch 
flight director signal). Similar behavior could be elicited In the roll 
and vertical axes by Increasing ULIM(2) and ULIM(3) threshold values. 

Finally, parameters ULIM(I) could be scheduled to vary with range to 
go from the helipad or be made adaptive to system performance levels (e.g., 
Instantaneous groundspeed, localizer and glldeslope errors) to reflect 
plausible human pilot decision and control strategies. 

3. Strategy Option ISTRAT = 1 

This option allows for closed loop piloted control of only one axis 
at a time. At any given Instant the pilot must decide to control either 
of the pitch, roll or vertical axes or none of the three based upon some 
Internal criterion for evaluating the Individual axis urgency for atten- 
tion and control. Here the pilot :s assumed to compare three decision 
functions a(l) * pitch flight dir.:tor signal, a(2) * roll flight direc- 
tor signal and .(3) * collective flight director signal and choose the 
axis with the largest decision function magnitude. Simulations with vari- 
ous reasonable values for a(I), I = 1,2,3 (e.g., such that 1" pitch = 1", 
roll =1" collective flight director signals which is the same as saying 
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Figure 28. NCC Program Simulation Results: CSS/FD/WIND/DH * 0 ft 
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Figure 28. NCC Program Simulation Results: CSS/FD/WIND/DH ■ 0 ft 
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Figure 28. NCC Program Simulation Results: CSS/FD/WIND/DH * 0 ft 
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ULIM(l) 

ULIM(2) 

l'LIM(3) 
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ISTRAT • 2» I CARLO - 10, 


0.134 In (■ 2 deg) 
0.03 In (■ 2 deg) 
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15 # pitch ■ 66* roll ■ 20 ft/s collective flight director signal) were 
Investigated. Results show that the overall system In this mode goes 
unstable In the first 30 seconds of flight. However, It would be pre- 
mature to rule out this operating strategy on the basis of this limited 
evaluation. Perhaps the a's need to be made adaptive to the nature of 
the specific phase of the approach to hover task. 


Determination of Pilot Acceptance 

The results above demonstrate the feasibility of using the parameters 
of a suitable pilot model as measures or metrics for relative task diffi- 
culty and hence pilot workload. Thus for the option 0 strategy, the maxi- 
mum permissible values for the control sampling Intervals INPTCH, INROLL, 
INCOLL in the pitch, roll and collective axis while maintaining acceptable 
performance are shown to be 7, 15 ana 20 seconds, respectively. Analysis 
of this type must be carried out for the other competing configurations 
to determine the maximum permissible values for INPTCH, INROLL and i XL 
for the same fixed level of acceptable system performance. These three 
numbers may then be compared Individually or In some aggregate form (such 
as the weighted sum of the sampling frequencies) to arrive at some measure 
of pilot preference or opinion rating. 

More work needs to be done on the problems associated with aggregating 
components of a workload vector metric Into a scalar workload measure. 
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CONCLUDING REMARKS 


A systematic approach towards the development of an analytical 
methodology for determining the level of pilot acceptance for a given 
system configuration Is formulated. A general computer program (NCC) 
Incorporating a nonlinear model for the UH-1H helicopter and a generic, 
homomophlc. dynamic model for the human pilots Information acquisi- 
tion, processing and control strategy was developed. 

The study points out that the weakest elements In the successful 
Implementation of the overall approach are the particular structure of 
the generic pilot model and the nature of the off-line simulation scheme 
for evaluating pilot-aircraft system performance. Specifically, the 
results of this investigation show that a generic pilot model capable 
of describing the human Information processing, decision making and con- 
trol strategies peculiar to the helicopter decelerating approach to a 
hover task Is necessary. In addition, an off-line pilot-aircraft 
simulation, Incorporating realistic models for the aircraft and the 
human pilot Is essential for a systematic analysis of system perfor- 
mance and pilot workload factors and their Impact on overall level of 
pilot acceptance. 

The UH-1H approach to a hover task is dominated by flight segments 
such as glide slope capture and deceleration to a hover that require ex- 
plicit deterministic pilot control responses to discrete step changes in 
the commanded flight conditions. As a result, a pilot model with the 
ability to describe predominantly transient response behavior as opposed 
to tracking perturbations about a nominal must be considered. Thus 
approaches which hypothesize the partioning of the pilot response into a 
nominal and a perturbation control strategy are Inappropriate. 

Finally, the study recommends the inclusion of reasonably accurate 
aircraft nonlinear equations of motion in the off-line simulation, and 
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the use of Monte Carlo methodology for obtaining approach path dispersions 
(1.e. ( ensemble means and standard deviations) of key system variables 
with range (-dlstance-to-go from the helipad) as the Independent variable. 
Linear, time varying covariance propagation methods are not recommended 
for tasks, such as decelerated approach to a hover, where the dominant 
response Is the deterministic nominal component, which must be obtained 
through a fully nonlinear simulation (note that the nominal response 
history cannot be obtained by sVng ng together a sequence of dynamic 
trim outputs). 

In summary, a structured approach towards evaluating the tradeoff 
between system performance and pilot workload and their relationship to 
overall pilot acceptance as expressed by a pilot opinion rating has been 
outlined. Further analytical and experimental efforts aimed towards the 
Identification of these reciprocal relationships and the development and 
validation of aggregated, scalar, performance and workload measures from 
truly multi attribute vector quantities are recommended. 
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APPENDIX: Outline of NCC Program 
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Figure Al. Flowchart for NCC 
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Figure A3. Listing of Subroutine F2 
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Figure A3. (Continued) 
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